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Abstract 

The lack of polar symmetry in two-bladed rotors 
leads to equations of motion with periodic coeffi- 
cients in axial flight, which is contrary to three or 
more bladed rotors that result in constant coefficient 
equations. With periodic coefficients, the analysis 
becomes involved, as a result very few studies have 
been directed towards the analysis of two-bladed ro- 
tors. In this paper, the aeroelastic stability of two- 
bladed proprotor/py Ion/ wing combinations is exam- 
ined in high speed axial flight. Several parametric 
studies are carried out to illustrate the special nature 
of two-bladed proprotors and to better understand 
the mechanism of whirl-flutter in such rotors. The 
wing beam bending mode for two-bladed rotors is 
found to be stable over the range of parameters ex- 
amined, a behaviour very different from three- bladed 
rotors. Also, the wing torsion mode exhibits a new 
type of instability similar to a wing torsional diver- 
gence occuring at 1/rev frequency. This type of be- 
haviour is not seen in three and more bladed ro- 
tors. The interaction between wing chord wise bend- 
ing and torsion modes is found to be much greater 
in the case of two-bladed rotors and, over the range 
of parameters considered, these two modes govern 
the stability of the system. 
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1 Introduction 

Aeroelastic st ability of a tiltrotor aircraft in forward 
flight has hei n the subject of considerable research 
since the 1960^s. In high speed axial flight, tiltro- 
tors are subject to an instability known as propro- 
tor whirl flutter. Due to high inflow through the 
rotor in this flight mode, large rotor forces depen- 
dent on the in-plane motion are generated. These 
forces are larger in tiltrotor aircraft than in conven- 
tional propeller driven aircraft due to large blade flap 
and lag motions in case of tiltrotors. With increas- 
ing airspeed, the rotor forces and the pylon/wing 
motion interact with each other to a point where 
the rotor forces become destabilizing and the ro- 
tor/pylon/wirig system becomes unstable. 

Most of the rtisearch on proprotor whirl flutter has 
focused on three or more blades. Several inves- 
tigations of this phenomenon were conducted us- 
ing pylon pivot models. Reed^^^ presented a re- 
view of these early efforts and highlighted some 
of the key issues associated with this instability. 
Johnson® developed a comprehensive mathemati- 
cal model of a tiltrotor including a modal repre- 
sentation of the wing, and for high speed forward 
flight. Analytical results obtained with this analy- 
sis showed good correlation with full scale proprotor 
tests £Lnd the analysis was later extended to include 
elastic blades and helicopter and conversion modes 
of operation. Kvaternik and Kohn^^^ carried out 
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an experimental parametric investigation of whirl- 
flutter for a proprotor mounted on a rigid pylon with 
flexibility in pitch and yaw. Results of this study 
showed that proprotor whirl flutter can be predicted 
with linear stability analyses using two-dimensional, 
quasi-steady aerodynamics for the blade loading. 
Nixon conducted parametric studies for aeroe- 
lastic stability of a tiltrotor in forward flight using a 
finite element based comprehensive tiltrotor aeroe- 
lastic analysis. The effects of several design parame- 
ters such as rotor flap and lag frequencies, pitch-flap 
coupling, wing stiffness and sweep were examined 
to determine their effect on aeroelastic stability in 
the high speed axial flight mode. The rotor was 
assumed to be in a state of autorotation, i.e., ro- 
tor torque was not transferred to the wing and the 
wing vertical bending slope at the tip did not af- 
fect rotor inplane motion. Johnson showed that 
this assumption gives conservative estimates of whirl 
flutter stability. Lag frequency tuning was found to 
be a practical method for increasing flutter speed in 
axial flight and forward sweep of the wing was found 
to be destabilizing (Nixon The above models 
rely on experimental data to model the lag frequency 
variation with blade collective pitch and the distri- 
bution of blade flexibility inboard and outboard of 
the pitch bearing. To reduce the dependance on ex- 
perimental data, Hathaway and Gandhi have re- 
cently made some refinements to the model which 
allow first principles based modeling of the lag fre- 
quency variation and blade flexibility distribution. 

In recent years numerous approaches for the im- 
provement of tiltrotor whirl mode boundaries for 
three-bladed tiltrotors have been investigated. Srini- 
vas and Chopra^^, Popelka^®^ et al., Corso^ et al. 
and, Barkai and Rand^^^ investigated the use of 
aeroelastic tailoring of composite wing and rotor 
blades to enhance whirl mode stability of a tiltrotor 
in forward flight. Nixon^^^^ et al. presented a review 
of the experimental and analytical efforts for sta- 
bility augmentation of tiltrotors through aeroelastic 
tailoring of the wing and rotor blades. Srinivas^^^ 
et al. carried out studies on improving whirl mode 
stability of tiltrotors by using advanced geometry 
blades with tip sweep and tip anhedral/dihedral. 

Acree^^^^ et al. investigated the effect of the chord- 
wise positions of the rotor blade aerodynamic center 
and center of gravity on stability. Rearward offeets 
of the aerodynamic center with respect to the blade 
elastic axis and pitch axis were found to give large in- 


creases in the stability boundary, with the effect be- 
ing most dominant for offeets at the outboard part 
of the blade. Acree^^^^ also analyzed the effect of 
blade tip sweep on whirl mode stability of the V-22 
tiltrotor with a tip mass extended on a boom for- 
ward of the leading edge to compensate for the aft 
shift of the tip center of gravity due to sweep. 

All of the above studies were carried out for three- 
bladed proprotors, although Johnson presented 
some results for a two-bladed propeller as well. Un- 
like a rotor with three or more blades, a two-bladed 
rotor lacks polar symmetry, which leads to propro- 
tor equations of motion containing periodic coeffi- 
cients. Thus, the stability analysis in the case of a 
two-bladed rotor is more involved and as a result, 
a lot less attention has been devoted to two-bladed 
configurations. In June 1998, whirl flutter was en- 
countered during a flight test of the Pathfinder air- 
craft, fitted with several two-bladed propellers to be 
used on the Helios aircraft. These aircraft, built un- 
der the Environmental Research Aircraft and Sen- 
sor Technology (ERAST) program, have lightweight 
and highly flexible wing structures and soft-mounted 
propellers. The simplified analysis for two-bladed 
rotors, when applied to the Pathfinder configura- 
tion with th»? Helios propellers, was unable to predict 
the flutter encountered in 1998. Also, it is unclear 
whether design guidelines established for preventing 
whirl flutter instability for three or more blades can 
be applied to the two-bladed case also. The peri- 
odicity of the equations of motion for a two-bladed 
proprotor may introduce new phenomena that re- 
quire careful examination and scrutiny. Hence a sys- 
tematic study of whirl flutter of two-bladed propro- 
tors is carried out to determine the effects of various 
design parameters on their stability. 

2 Analytical Model 

Figure 1(a) .shows the coordinate systems used for 
the rotor and wing. The fixed frame rotor forces, 
which are tr.ansferred to the wing, are also shown. 
Figure 1(b) shows the geometry of the wing/pylon 
system. Flap and lag motions of the rotor blades are 
considered and, the blades are assumed to be rigid. 

2.1 Rotor and Hub Structural Model 

The equations of motion for a single rigid blade 
in the rotating frame with flap zmd lag degrees of 
freedom have been presented by Johnson These 
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(a) Coordinate Systems for Rotor and Wing 



(b) Geometry of the Wing and Pylon 


Figure 1: Coordinate Systems for Rotor /Wing and 
Wing Geometry 

equations include the effects of pylon translation and 
rotation on the flap and lag motions^ 

+ 2ay)sinil>m] + S'pZp = 7 " - " - ( 1 ) 

("’) + + SI [(i*p + hay)sinrl)m - (Vp 

- haj:)cosrf>m} - = 7 — — (2) 

where the flap motion, lag 

motion and the azimuthal position of the blade, 
and are the non-dimensional flap and lag fre- 
quencies, respectively, and h is the distance between 
the pylon pivot point and the hub. Xp, j/p and Zp are 
pylon translational degrees of freedom, as shown in 


Figure 1, while a*, ay and a* are pylon rotational 
degrees of freedom. and are respectively 

the aerodynamic flap and lag moments of the 
blade. The above equations can be written in matrix 
form as, 

** ^ * Mjt 

\^^]Pr + + \^^p\p — (^) 

(/c‘l?H + (^C*1C« + l^cVl*^ = (4) 

where Pr and Cr *u:e vectors of the flap and lag 
motions in the rotating frame while p contains the 
pylon degrees of freedom. For N blades, 

0R = [P^^^ ... 

C« = [cw -- 

p~-[xp Vp Zp Ox Qy 

Other matrices in Equations 3 and 4 are defined in 
the Appendix. For coupling with the pylon/wing 
system, the above equations must be transformed to 
the fixed frame. Also, the pylon degrees of freedom 
must be replaced by the wing-tip degrees of freedom. 
This is accomplished using the following transforma- 
tions, 

Pr = [Trf\Pf 6e = P = [T\]Wtxp 

where pp and are the fixed frame rotor degrees of 
freedom, while Wxip is a vector containing the finite 
element degrees of freedom at the wing- tip. 
is an iV X A' matrix relating the fixed frame rotor 
degrees of freedom to the rotating frame degrees of 
freedom, while |Ta] is a transformation matrix that 
depends on the wing sweep A and relates the pylon 
degrees of freedom to the wing-tip degrees of free- 
dom. Ppy Cfi [^rf] and [Ta] are as defined in the 
Appendix. 

For a two-bladed rotor, the collective flap (/^), dif- 
ferential flap (Pi)y collective lag (Co) and diflerential 
lag (Cl) motions comprise the fixed frame rotor de- 
grees of freedom. In case of a three-bladed rotor, the 
fixed frame motions are collective flap (A)), cyclic 
flap {pic and A,), collective lag (Co) and cyclic lag 
(Cic and Cij)- 

Note that, at this stage, the number of blades is 
assumed to be arbitrary so that the same analysis 
can be applietl to both two and three-bladed rotors. 
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(a) Velocity and Force Components on a Typical 
Blade Section 



Figure 2: Blade Section Velocities and Wing Ele- 
ment 


Sup = r(/l - ayCOsi)m + 


= r6upg + Jup^ 


( 6 ) 


6ur — sinV'm - Xpcosrpm + h{axSin^l>m - 


COSrl>„ )- - CtyCOS^rn) 


(7) 


In Equations 5,6 and 7, Suxj^t 

Sur are independent of the radial position r and 
hence they may be factored out of the integrands of 
the aerodynamic force integrals. The aerodynamic 
moments on a blade can then be written in terms of 
the aerodynamic coefficients (M^, etc.) and the 

perturbation velocities (see Johnson^^). For exam- 
ple, the perturbation flap moment may be written 
as, 




-h MffSe (8) 

where SO is the perturbation in the collective angle. 
The perturbation velocities, Suta^ ^^b 
functions of 0rXR} ^Tip and ipm for each blade, 
and thus the right hand side terms of Equations 3 
and 4 can be expressed in terms of the system de- 
grees of freedom and moved over to the left hand 
side of the etjuations. 


2.2 Rotor Aerodynamics 

Figure 2(a) shows the velocity components at a blade 
section defined with respect to the hub plane. Aero- 
dynamic forces are based on a linear strip theory 
and the rotor is assumed to be in a state of autoro- 
tation. The trim procedure is based on changing 
the collective pitch until the thrust becomes zero. 
Perturbations in the blade section velocities, for a 
forward flight speed V, can be written as, 

Sut - r(az - C) + {-XpSiniprn + VpCOSlprn - 
* V 

COStprn siniprn ) + {o^x COSlprn~^ay sin^l/}^)) 


= tSuta + (5) 


2.3 Wing Structural Model 

The wing structural model is based on a finite el- 
ement formulation of an Euler-Bernoulli beam un- 
dergoing beamwise bending, chordwise bending and 
torsion, as shown in Figure 2(b). The axial degree of 
freedom is not considered since the wing is assumed 
to be rigid in axial direction. 

2.4 Wing Aerodynamics 

Wing aerodynamics is based on a quasi-steady 
approximation and the aerodynamic stiffness and 
damping matrices are determined and added to 
the corresponding structural matrices. Wing sweep 
(Figure 1(b)) is taken into consideration in formu- 
lating the wing aerodynamic matrices. 
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2.5 Coupling of Rotor/Wing System 

Coupling of the rotor/ wing system essentially con- 
sists of two parts: (1) incorporating the rotor aero- 
dynamic and inertial forces in the fixed frame into 
the wing model, and (2) incorporating the wing-tip 
degrees of freedom into the rotor equations of mo- 
tion (see Sec. 2.1 and 2.2), this leads to a system of 
equations with periodic coefficients. In state space 
form, this system can be written as, 

y = A{^)y (9) 

♦ 

where ( ) denotes a derivative with respect to the 
azimuth, — A{xIj) is a periodic matrix with 

period 27 T and y is a column vector of the 2 m system 
states. 

2.6 Stability Analysis 

The stability of the system defined by Equation 9 is 
examined using Floquet analysis. The Floquet tran- 
sition matrix $(27r), relates the states of the system 
at ^ = 0 to the states at V* = 27 t, 

y(27T) = 4>(27T)y(0) (10) 

Each column of the Floquet transition matrix 
(FTM) is determined by integrating Equation 9 with 
a unit initial condition for one state and zero ini- 
tial conditions for other states. Thus 2m integra- 
tions need to be performed to obtain the FTM, 
where 2m refers to the number of states. An effi- 
cient method, developed by FViedmann et.al.t^^^, for 
computing the FTM in a single integration pass is 
used in the present study. Eigenvalues of the FTM 
are called Floquet multipliers and are denoted by 
di (i=l,2...2m). The system is unstable if any of 
the Floquet multipliers have an absolute value larger 
than unity (see Johnson^^^l, Nayfeh and Mook^^^). 
Floquet multipliers are related to the Floquet expo- 
nents (Aj) as, 

A, = ^ln(ei) (11) 

Z7T 

Real and imaginary parts of the Floquet exponents 
are related to the conventional damping (Q) and fre- 
quency (ctfpi) respectively. 

O = = (12) 

Ljpi = Jm{Xi) = (13) 


^Co 

Rotating Blade 
Freq. (per rev) 
In- vacuo 
1.02 
1.85 

see Fig. 3 
0 


Wing FVeq. (per rev) 

Beam 

0-42 

(yhord 

0.70 

Torsion 

1.3 


Table 1: FVequencies of the Baseline System 



Figure 3: Blade Lag FVequency as a Function of Ve- 
locity 

A subscript ’p* for the frequency indicates that only 
a principal value of the frequency can be determined 
due to the fact that the logarithm of a complex num- 
ber has many branches, giving values for Aj that 
differ in frequency by integer multiples of the rota- 
tional speed Q. Compared to a constant coefficient 
system, the periodic system is also described by nor- 
mal modes (i.e, eigenvectors) and the roots A^, but 
for a periodic system the eigenvectors are periodic 
functions rather than constants (see Johnson ^^^^). 

The analytiad model was validated by comparing 
the frequency and damping variations with flight 
speed, with predicted results obtained by Nixon^^^ 
for a three b laded rotor. Also, the whirl flutter 
stability boundaries for a two-bladed propeller were 
compared with the results reported by Johnson^^J. 
In both cases, the present results matched with those 
of Nixon and Johnson. 
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(a) Frequency 



(b) Damping 


Figure 4: Frequency and Damping as a Function of 
Velocity for Baseline System (Two-Bladed Rotor) 

3 Baseline System 

The present model is used to obtain the frequency 
and damping characteristics of two and three-bladed 
rotor/wing systems. The baseline rotor/wing system 
has the same properties as the baseline three-bladed 
rotor used by Nixon No wing sweep is considered 
in the present study. Table 1 shows the important 
in- vacuo frequencies of this system. The collective 
flap mode has a higher frequency because the blades 
act as cantilevers in this mode due to a gimballed 
hub. The in- vacuo lag frequency varies with collec- 
tive pitch (and therefore flight speed) due to the dis- 
tribution of lag flexibility inboard and outboard of 



(a) FVequency 



(b) Damping 


Figure 5: FV^iquency and Damping as a Function of 
Velocity for Baseline System (Three-Bladed Rotor) 


the pitch bearing (Figure 3). The rotor has a radius 
R=12.5 ft, Lock number 7=3.83, solidity a=0.089, 
pitch- flap coupling =-0.268 and rotational speed 

n=48 rad/s. The wing/pylon system has a span of 
16.65 ft, a chord of 5.16 ft and a pylon mast height 
of 4.28 ft. I’hree finite elements of length 4.55 ft 
each are used to model the wing, with one pylon el- 
ement of length 3.0 ft at the tip. The wing has a 
beamwise stiffness Eh=3.l3 x 10^ Ib-ft^, chordwise 
stiffriess EIc~SAS x 10^ and torsional stiffness 
GJ=1.62 X 10^ ll>ft2. 

Figure 4 shows the frequency and damping charac- 
teristics of the baseline two-bladed rotor. It must be 
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noted that the modes are highly coupled and hence 
their identification is based on the physical frequency 
and the behavior of periodic eigenvectors. In Figure 
4, *b’, ’c’ and ’t’ refer to modes that are primarily 
wing beam, chord and torsion modes respectively, 
while ’/3b, and ’Ci’ are the rotor collective flap, 
differential flap and differential lag modes. The same 
nomenclature is used in all subsequent figures. 

Shown in Figure 4 are the frequency and damping 
characteristics computed using two methods: (1) us- 
ing a full 48 state system with 3 wing elements and 
1 pylon element representing the wing, and (2) us- 
ing an 18 state system obtained by applying static 
condensation to the wing equations and retaining 
only the five wing-tip degrees of freedom, prior to 
coupling with the rotor equations. The reduced sys- 
tem captures all the essential modes quite well and 
is computationally much faster than the full system. 
Further results in this paper are based on the re- 
duced order system . 

Figure 5 shows the frequency and damping char- 
acteristics of a three- bladed rotor having the same 
non-dimensional properties as the two-bladed rotor. 
As compared to the two-bladed rotor, two addi- 
tional degrees of freedom are introduced in the three- 
bladed rotor due to the flap and lag motions of the 
third blade. Thus for a three- bladed rotor, there 
are five important rotor modes (/^, 4* 1, — 1, 

4- 1, ^ — 1) compared to three rotor modes(^,/3i 
and Cl) for the two-bladed case. In Figure 5, — 1 

and C “ 1 refer to the low frequency flap and lag 
modes respectively, while ^4-1 and C “h 1 are the 
high frequency flap and lag modes and j3o Is f be col- 
lective flap mode. The collective lag mode Co Is not 
shown in Figures 4 and 5 since it has zero frequency 
due to the assumption that the rotor is in a state of 
autorotation. 

A key difference between two and three- bladed ro- 
tor behaviour can be seen in the wing beam mode 
damping for the two cases. For the baseline three- 
bladed rotor, the wing beam mode becomes unstable 
for V/VlR > 1.19 as seen from the beam mode damj>- 
ing in Figure 5(b). However, it is evident from the 
beam mode damping for a two-bladed rotor (Figure 
4(b)) that the wing beam mode remains stable for 
all flight speeds considered. An explanation of this 
behaviour is presented in Section 4.1. 

For a three- bladed rotor, there is an energy transfer 


between the wing chord mode c, and the low fre- 
quency lag mode, C “ 1, at the point where the C “ 1 
frequency crosses the chord mode frequency (Figure 
5(b), V/UR = 0.2). The chord mode damping re- 
duces slightly at this point since the damping in the 
C — 1 mode is lower. A similar interaction is seen 
when the C “ 1 niode frequency crosses the beam 
mode frequency at V/Q,R = 0.4 for the three-bladed 
rotor (Figure 5), where the beam mode damping 
shows a peak since, at this point, the damping in 
the C ” 1 niode is higher than that in the beam 
mode. Figure 4 presents evidence of similar energy 
transfer between the differential lag mode Ci and 
the wing beam and chord modes for a two-bladed 
rotor. This occurs at the points where the princi- 
pal value of the Ci mode frequency crosses the wing 
chord and besun mode frequencies at V/VtR = 0.22 
and 0.48, respectively (Figure 4(a)). As compared 
to a three-bladed rotor, the energy transfer between 
the rotor lag and wing modes is smaller in the case of 
a two-bladed rotor as indicated by the smaller peak 
in the beam mode damping for the two-bladed ro- 
tor (Figure 4(b)). Also, the valley in the chord mode 
damping at V/hR — 0.22 (Figure 4(b)) is very small 
for the two-bladed rotor. Unlike a three-bladed ro- 
tor the frequencies need not actually cross in case of 
a two-bladed rotor (i.e., only their principal values 
cross) since, for a two-bladed rotor, Floquet mul- 
tipliers govern the interaction between the modes. 
For the two-bladed rotor the energy transfer occurs 
when both the real and imaginary parts of the Flo- 
quet multipliers of two modes approach each other. 

In the case of a two-bladed rotor, for V/HR > 1.1, 
the Floquet multiplier for wing torsion becomes real, 
which causes a split in the damping curve for this 
mode, one part increasing and the other decreasing 
to become unstable at V/HR > 1.25 (Figure 4). This 
is a unique instability for a two-bladed rotor, similar 
to the 1/rev divergence of a two-bladed propeller 
(Johnson f^). The differential lag mode also shows 
similar behavior at a high flight speed (Figure 4, 
V/QR > 1.9). but remains stable. 

4 Parametric Studies 

Several parametric studies are carried out to deter- 
mine the effect of key natural frequencies of the ro- 
tor and wing on the critical flight speed and to ex- 
amine the interactions between various modes for a 
two-bladed rotor. The critical flight speed Vf is the 
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Figure 6: Effect of Flap Frequency (i/p) on Critical 
Flight Speed 

speed at which the damping for a particular mode 
becomes zero. For comparison, results of these para- 
metric studies are plotted along with the results for a 
three bladed rotor having the same non-dimensional 
properties as the two-bladed rotor. 

4.1 Rotor Frequencies 

The rotor frequencies, up and are independently 
varied from their baseline values by factors varying 
between 1.0 and 2.2. It must be noted that extreme 
variations of in-vacuo flap and lag frequencies may 
not be practically possible, but they are of theoret- 
ical importance to establish the trends within the 
design range. Also, since the baseline lag frequency 
is a function of velocity as shown in Figure 3, the 
factor for shifts the entire curve by a factor be- 
tween 1.0 and 2.2. Since up and u^^ are uncoupled, 
in-vacuo flap and lag frequencies of the blades, an 
increase in these frequencies must primarily effect 
the pi or Cl frequencies which are the differential 
flap and lag frequencies of a coupled system with 
aerodynamics. 

Figures 6 and 7 show the results of varying up and 
u^ by the indicated factors from their baseline val- 
ues. Dotted lines show the results for a three-bladed 
rotor having the same non-dimensional properties 
as the two-bladed rotor. Note that, in these fig- 
ures, only the most critical modes for whirl stability 
are shown. Figures 8 and 9, respectively show the 
frequency and damping characteristics for two and 
three-bladed rotors when up is increased by a factor 
of 1.1 from its baseline value. Frequency and damp- 
ing characteristics, with up increased to 1.5 times its 



Figure 7: Effect of Lag Frequency (i/^) on Critical 
Flight Speed 

baseline value, are shown in Figures 10 and 11 for 
two and thri^e-bladed rotors respectively. The prin- 
cipal observiitions (italicized) are enumerated below 
along with their explanations (non-italicized). 

1. For a two-bladed rotor, the wing beam mode b does 
not become ^mstable for the covered range of up and 
U(^. This is a key difference between two and three- 
bladed rotors. It is evident from Figures 6 and 7 
that the wing beam mode is quite critical for three- 
bladed rotors, especially for low flap frequency, but 
it is stable fcir the two-bladed case. 

For a three-])laded rotor, the wing beam mode be- 
comes unstable due to its interaction with the p—1 
and torsion modes (see Nixon^^^). As shown in Fig- 
ures 5 and 9, the C ~ 1 mode frequency acts as a 
barrier between the p ^ \ and beam mode frequen- 
cies at low flight speeds, but when the ^ — 1 fre- 
quency crosses above the C “ I frequency it drives 
the beam m(»de into an instability. As up factor is 
increased bevond 1.4, the beam mode for a three- 
bladed rotor does not become unstable as shown in 
Figure 6. N;xon^^^ showed that this is because of 
increased interaction between the rotor flap and lag 
modes and due to an increase in the wing torsion fre- 
quency which moves it further away from the wing 
beam mode frequency. 

For a two-bladed rotor the differential flap {p\) fre- 
quency does not interact strongly with the wing 
beam mode frequency for the covered range of up 
as shown in Figures 4, 8 and 10. Thus, the wing 
beam mode remains stable for a two-bladed rotor. 

2. As shown in Figure 6 the wing chord mode for 
a two-bladed ivtor initially becomes less stable as up 
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(a) FYequency 



(b) Damping 


Figure 8: Frequency and Damping Characteristics 
for Factor 1.1 (Two-Bladed Rotor) 

is increased and then becomes more stable. Chord 
mode damping for a two-bladed rotor is affected by 
the placement of the P \ , Ci and wing torsion frequen- 
cies relative to the chord mode frequency. The 
frequency exerts a destabilizing influence on wing 
chord mode damping, whereas wing torsion fre- 
quency has a stabilizing influence on it. For low i/p 
factors (i/p < 1.5), the mode frequency falls be- 
low the torsion mode frequency for V/QR ~ 0.5 — 1.0 
(see Figures 4 and 8). In this range, the chord mode 
damping reduces faster and the overall effect is a 
slight reduction in critical flight speed for up fac- 
tors 1.1 and 1.2. For up factors above 1.5, the Ci 



(a) FVequency 



(b) Damping 


Figure 9: FVequency and Damping Characteristics 
for up Factor 1.1 (Three-Bladed Rotor) 

mode frequency does not cross below the torsion 
mode frequency and also the mode frequency 
moves further away from the chord mode frequency 
for V/flR < 1.1 as shown in Figure 10. This stabi- 
lizes the wing chord mode at low flight speeds, but 
the increase in its critical flight speed is not signif- 
icantly large since the 0i mode frequency interacts 
very strongly with the chord mode at higher flight 
speeds (Figure 10(a)). 

In case of a three-bladed rotor the wing chord mode 
gains greater stability as is increased, as shown 
in Figure 6 since, for a three-bladed rotor, the in- 
stability in the wing chord mode is caused by its 
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(a) Ftequency 



(b) Damping 


Figure 10: FVequency and Damping Characteristics 
for Factor 1.5 (Two-Bladed Rotor) 

interaction with the p — 1 mode (Nixon As i/^ 
is increased, the p—1 mode frequency crosses above 
the C “ 1 mode frequency at a higher flight speed 
(see Figures 5, 9 and 11). This leads to an increase 
in the critical flight speed for the wing chord mode 
in case of a three-bladed rotor. 

FYom the above discussion it is clear that the wing 
chord mode is driven into an instability due to its 
interaction with the rotor flapping mode, for both 
two and three-bladed rotors and, in both cases, the 
placement of the rotor lag frequency is important. 
However, in case of a two-bladed rotor there is a sig- 
nificantly greater interaction between the wing chord 



(a) Frequency 



(b) Damping 


Figure 11: Frequency and Damping Characteristics 
for Factor 1.5 (Three-Bladed Rotor) 

and torsion modes. Further evidence of this interac- 
tion is presented in Section 4.2. 

3. As shown in Figure 6, the wing chord mode for a 
two-bladed rotor does not become unstable when up 
factor is greater than 2My but an instability develops 
in the mode at this point For a two-bladed rotor, 
as the p\ frequency increases, it has a lesser inter- 
action with the chord mode and greater interaction 
with the Cl niode. Consequently, Ci mode damping 
reduces sharply as shown in Figures 4, 8 and 10. 

For a three-bladed rotor also, there is a significant re- 
duction in the damping for the ( — 1 and C + 1 modes 
due to greatei interactions with the rotor flap modes 
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(a) Frequency 



(b) Damping 


Figure 12: FVequency and Damping Characteristics 
for Factor 1.1 (Two-Bladed Rotor) 

(see Figures 5, 9 and 11), as vp is increased. Such 
interaction leads to an instability in the highest fre- 
quency rotor mode (Nixon ^^^). As shown in Figure 
6, for factors between 1.3 and 1.95, the C+1 mode 
is unstable in case of a three-bladed rotor. When 
factor increases above 1.95, the + 1 mode becomes 
the rotor mode with the highest frequency and thus 
the instability shifts to the ^ -h 1 mode (Figure 6). 

For a two-bladed rotor, the instability in the Ci mode 
occurs only when the factor is increased beyond 
2,0, in contrast to the three-bladed case where the 
(-hi mode becomes unstable for factors greater 



(a) FVequency 



(b) Damping 


Figure 13: frequency and Damping Characteristics 
for Factor 1.1 (Three-Bladed Rotor) 

than 1.3 (Figure 6). This is because for a two- 
bladed rotor, at high values of the 01 mode fre- 
quency has a close interaction with the (i mode fre- 
quency only at low flight speeds (see Figure 10). At 
higher flight speeds, the 0\ mode frequency falls be- 
low 1/rev and interacts more with the wing chord 
mode. This stabilizes the (i mode at higher speeds 
as shown in Figure 10(b). For a three-bladed ro- 
tor, the 0 + \ mode frequency is close to the ( + 1 
mode frequency even for the baseline case (Figure 
5). Thus, only a small increase in yp is needed to 
increase the interaction between the C “b 1 + 1 

modes and thus drive the C + 1 mode unstable in 
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(a) FVequency 



(b) Damping 


Figure 14: FVequency and Damping Characteristics 
for Factor L5 (Two-Bladed Rotor) 

case of a three- bladed rotor, whereas a much larger 
increase in is required for the 0i mode frequency 
of the two-bladed rotor to cause an instability in the 
(i mode. 

4. As shown in Figure 6 the wing torsion mode 
becomes more stable as up factor increases and be- 
comes completely stable for up factors above 1.6. For 
a two-bladed rotor, as the 0i frequency increases, 
it causes an increase in the wing torsion frequency 
as well (see Figures 4, 8, 10). This increases the 
flight speed at which the wing torsion mode diverges 
(i.e., speed at which the Floquet multiplier for tor- 
sion becomes real) and thus there is an increase in 



(a) FVequency 



(b) Damping 


Figure 15: Frequency and Damping Characteristics 
for u^ Factoi 1.5 (Three- Bladed Rotor) 

the critical flight speed. Also, overall increase in 
the 01 frequency causes it to rise above the chord 
mode frequency for V/QR > 1.2 (see Figure 10), 
and thus interact more with the torsion mode. This 
causes an increase in damping for the torsion mode 
at higher flight speeds, ultimately making this mode 
completely stable. 

For a three-bladed rotor also, the wing torsion mode 
becomes more stable as up is increased as seen from 
the torsion mode damping in Figures 5, 9 and 11. 
This is because of an increase in torsion mode fre- 
quency which moves it further away from the 0 
mode frequency. 
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5. For a two-bladed rotor, the wing chord mode ini- 
tially becomes less stable as factor is increased to 
1.2, then becomes more stable with factor and re- 
mains relatively unaffected for factors above 1.8 
(Figure 7). There are three effects at work here. 
Firstly, an increase in the Ci frequency causes an in- 
crease in the flight speed where its principal value 
crosses the chord mode frequency. Thus the cross- 
over now takes place at a point where the Ci mode 
damping is higher than the chord mode damping, 
resulting in a damping peak in the chord mode in- 
stead of a valley (Figure 12(b)). This peak shifts to 
higher flight speeds as factor increases. Secondly, 
an increase in frequency reduces the range of flight 
speeds over which it is below the torsion frequency 
(see Figures 4 and 12) . This is a stabilizing effect 
for the chord mode as discussed before. Thirdly, 
the portion of the !3\ frequency below the chord fre- 
quency increases so that there is greater interaction 
between the and chord modes, which is destabi- 
lizing (Figures 12(a) and 14(a)). For factors of 
1.1 and 1.2, the third effect is slightly greater than 
the first two leading to a slight decrease in the crit- 
ical flight speed- For higher factors, the first two 
factors cause an increase in the critical flight speed. 
For factors above 1.8, the Ci frequency does not 
have a great effect on any other mode and thus the 
critical flight speed for the chord mode becomes rel- 
atively constant under the opposing influences of the 
Pi and wing torsion modes. 

For a three- bladed rotor, the wing chord mode be- 
comes more stable as is increased, until about 
then remains relatively unaffected (Figure 7). 
For low values of the damping in the wing chord 
mode for a three-bladed rotor is strongly influenced 
by the location of the C - 1 frequency (Nixon 
As increases, there is an increase in the ^ — 1 fre- 
quency as seen in Figure 15. This increases the flight 
speed at which the C ~ 1 mode frequency crosses the 
chord mode frequency, leading to an increase in the 
critical flight speed for the chord mode. Beyond a 
factor of 1.7, the C — 1 frequency moves above the 
wing torsion frequency and thus the flutter in the 
chord mode becomes dominated by the ^ — 1 mode 
(Nixon ^^^). 

For both two and three-bladed rotors, at higher val- 
ues of i/^, the wing chord mode flutter becomes dom- 
inated by the flap mode. However, for low the 
chord mode for a two-bladed rotor behaves differ- 
ently due to the greater interaction between the wing 
torsion and chord modes for a two-bladed rotor. 


6. As shown in Figure 7, for a two-bladed rotor, 
the wing torsion mode stabilizes slightly as fac- 
tor increases to 1.7 and then destabilizes slightly 
for greater factors. For both two and three- 
bladed rotors, an increase in causes an increase 
in the torsion mode frequency which is most pro- 
nounced at low flight speeds (Figures 14(a) and 
15(a), V/VlR < 0.3). For a two-bladed rotor this 
leads to an increase in the flight speed at which the 
torsion mode diverges, i.e. the flight speed at which 
the torsion mode Floquet multiplier becomes real. 
This causes the split in the torsion mode damping 
to occur at a higher flight speed and thus there is an 
increase in t he critical flight speed. As the fac- 
tor increases beyond 1.7, the effect of frequency 
on the torsion frequency diminishes in the case of 
a two-bladeil rotor. This causes a slight decrease 
in the flight speed at which the torsion mode di- 
verges, leading to a slight decrease in the critical 
flight speed. 

4.2 Wing Frequencies 

Figures 16, 17 and 18 show the variation of criti- 
cal flight spt^ as the wing beam, chord and tor- 
sion stiffnesses (i.e. frequencies) are varied (indepen- 
dently) from their baseline values. Nixon showed 
that, for a three-bladed rotor, the stability of the 
wing beam mode depends on the proximity of the 
wing beam mode frequency to the wing torsion fre- 
quency As shown in Figure 16, when the wing 
beam stiffness is reduced for a three-bladed rotor, 
the beam mode becomes more stable because the 
beam mode frequency moves away from the torsion 
frequency Conversely, a reduction in the wing tor- 
sional stiffriess for a three-bladed rotor brings the 
torsion mode frequency closer to the beam mode 
frequency, leading to a reduction in beam mode sta- 
bility as shown in Figure 18. Nixon^^J also showed 
that varying the wing frequencies causes a mild ef- 
fect on wing chord mode stability for a three-bladed 
rotor because*, in this case, the wing chord mode is 
not highly coupled to the other two wing modes (see 
Figures 16, 17 and 18). 

For a two-bladed rotor, varying the wing beam fre- 
quency has a mild effect on the stability of the wing 
chord and torsion modes (Figure 16), since the beam 
mode does not have much interaction with any of 
the other modes. When the frequency of the wing 
beam mode increases, it increases the flight speed 
at which the Floquet multiplier for the wing torsion 
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Figure 16: Effect of Wing Beam Stiffness [Eh) on 
Critical Flight Speed 

mode becomes real. This leads to a small increase 
in the critical flight speed for the wing torsion mode 
of a two-bladed rotor. Also, greater interaction be- 
tween the beam and torsion modes slightly reduces 
the stabilizing influence that the torsion mode ex- 
erts over the chord mode in the case of a two-bladed 
rotor. This leads to a slight decrease in the critical 
flight speed for the chord mode as wing beamwise 
stiffness {Elb) is increased. 

The effects of wing chordwise stiffness {EIc) and 
wing torsional stiffness (GJ) on critical flight speed 
of a two-bladed system are more pronounced. Sev- 
eral observations regarding Figures 17 and 18 are 
enumerated below (italicized) along with their ex- 
planations (non-italicized). 

1. As shown in Figure 11 ^ the wing chord mode for a 
two-bladed rotor becomes more stable os its frequency 
is reduced. As discussed earlier, for a two-bladed ro- 
tor, wing chord mode damping depends on the place- 
ment of the pi , Cl and wing torsion frequencies with 
respect to each other and also with respect to the 
wing chord mode frequency. As the chord mode fre- 
quency reduces, the p\ and Ci frequencies move fur- 
ther away from the chord frequency for V/Q,R < 1. 
This increases the stability of the wing chord mode. 

2. For a two-bladed rotor, increasing the wing chord 
mode frequency from its baseline value leads to a 
slight reduction in its critical flight speed initially 
(up to EIc factor 1.2), and then a slight increase 
takes place (Figure 17). As the wing chord mode fre- 
quency increases it comes closer to the torsion and 
Pi frequencies for a two-bladed rotor. The net ef- 
fect is a slight decrease in stability initially, due to 



Figure 17: Effect of Wing Chordwise Stiffness {EIc) 
on Critical Flight Speed 



Figure 18: Effect of Wing Torsional Stiffness (GJ) 
on Critical Flight Speed 

the influence of the pi mode, and then an increase 
due to the interaction between the wing chord and 
torsion modiis. 

3. Increasing the wing chord mode frequency has a 
small stabilidng influence on the wing torsion mode 
of a two-bladed system (Figure 10). In this case, as 
EIc increases, greater interaction between the wing 
chord and torsion modes causes a slight reduction 
in the rate at which the torsion mode damping re- 
duces, beyond the point of divergence (i.e., the air- 
speed at which the Floquet multiplier for wing tor- 
sion becomes real). This causes a slight increase in 
the critical flight speed for the torsion mode. 

4. As shown in Figure 18, a reduction in the wing 
torsion frequency from its baseline value causes a 
large increase in chord mode critical flight speed for 
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(a) Frequency 



(b) Chord Mode Damping 


Figure 19: FVequency and Chord Mode Damping 
Characteristics for Various GJ Factors (Two-Bladed 
Rotor) 

a twO’bladed system. There are two effects at work 
here. Firstly, as the torsion mode frequency reduces, 
the Cl frequency is unable to cross below it. This 
causes a significant reduction in the rate at which 
the chord mode damping reduces with flight speed, 
due to greater interaction between the torsion and 
chord modes (see Figure 19). Secondly, a reduction 
in torsion frequency reduces the portion of the 
frequency that lies below the chord frequency (see 
Figure 19, V/QR >1). This leads to reduced inter- 
action between the Pi and chord modes, leading to 
an increase in chord mode damping at higher flight 



(a) FVequency 



(b) Chord Mode Damping 


Figure 20: FVequency and Chord Mode Damp- 

ing Characteristics for Various GJ Factors (Three- 
Bladed Rotor) 

speeds. The overall effect is a large increase in chord 
mode critical flight speed for the two-bladed system, 
as GJ is reduced. 

5. For a two-bladed roior, o reduction in the wing 
torsion frequency significantly reduces the stability 
of the torsion mode (Figure 18). This is primarily 
a result of the loss in stif&iess of the wing torsion 
mode, which reduces the airspeed at which the Flo- 
quet multiplier for wing torsion becomes real (i.e. 
the airspeed at which the wing torsion frequency re- 
duces to 1/rev). Also, as the torsion frequency re- 
duces, the interaction between (^i and torsion modes 
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reduces for low airspeeds {V/Q,R < 0.5). This causes 
a reduction in torsion mode damping at low speeds. 
Consequently, when the torsion mode frequency be> 
comes 1 /rev, it has a lower damping for low G J fac- 
tors. This causes a significant reduction in torsion 
mode stability. 

Similar to a three- bladed system, the wing beam 
mode damping for a two-bladed rotor also reduces 
when GJ is reduced. This is because the wing beam 
mode is coupled to the wing torsion mode through 
the wing chordwise mass offset. However, for a two- 
bladed rotor, the beam mode remains stable for most 
GJ factors, becoming unstable (at V/D.R = 1.24) 
only when GJ is reduced by half (not shown in Fig- 
ure 18). 

The most important result of the wing frequency 
study is the greater coupling between the wing chord 
and torsion modes for a two-bladed rotor in contrast 
to the three-bladed case. Figures 16, 17 and 18 show 
that wing chord and torsion modes are the governing 
modes for stability of a two-bladed proprotor. It is 
also evident that for a three-bladed rotor, changing 
the wing frequencies has little effect on chord mode 
stability. But for a two-bladed rotor, the effect of 
EIc factor and GJ factor on chord mode stability 
is quite pronounced. This indicates that there is 
greater interaction between chord and torsion modes 
for two-bladed rotors, as discussed in Section 4.1. A 
comparison of Wing chord mode damping for two and 
three-bladed rotors (see Figures 19(b) and 20(b)), as 
affected by the wing torsional frequency, confirms 
the increased coupling of wing chord and torsion 
modes for two-bladed systems. In Figure 19(a), for 
low flight speeds {V/Q,R < 1.0), only the wing tor- 
sion frequency is effected by the change in wing tor- 
sional stiffness. Thus, for the two-bladed system, the 
change in chord mode damping (Figure 19(b)) must 
be due to the interaction between the wing torsion 
and chord modes. In contrast, for the three-bladed 
system, the change in wing torsional frequency has 
only the slightest effect on chord mode damping 
(Figure 20), indicating a very weak coupling between 
wing chord and torsion modes. 

5 Conclusions 

The results of several parametric studies on the sta- 
bility of two-bladed proprotors are presented. These 
studies illustrate the special nature of two-bladed ro- 
tors as opposed to three or multi-bladed rotors. Un- 
like a three-bleided rotor, the wing beam mode for a 
two-bladed rotor does not interact much with any of 


the other modes £ind thus remains stable. The wing 
torsion mode exhibits a unique instability similar to 
a wing torsional divergence at 1/rev frequency. Also, 
the interaction between the wing chord and torsion 
modes is much more pronounced in the case of two- 
bladed rotors. It is shown that the stability of the 
wing chord mode is governed by a complex interac- 
tion between the rotor differential flap, differential 
lag and the wing torsion and chord modes. Also, the 
wing chord and torsion modes are the two modes 
that govern the stability of two-bladed proprotors 
over the range of parameters considered. 
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Appendix 

For an even number of blades, 

= [/So 01c 0U 0Mc 0Ma 0N/2]^ 

( 1 ) 

6^ = [Co Clc (is (Mc (Ms Ca^/2]^ (2) 

For an odd number of blades. 


0F = 

[A) 01e 0\a • • • 0Mc 0Ms\ 

(3) 

(f = 

[Co Clc (la • • • (mc (Ms\ 

(4) 

where, 
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where, 

Sm = sinipm and Cm = cosipm 
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The column [—1 1 ... (—1)^] must be ap- 

pended to [Thf] if N is even. 

The wing-tip degrees of freedom are, 

Wnp = [wt Vi v't 4>tf (14) 

Pylon degrees of freedom are given by, 

p= [TA]WTip ( 15 ) 
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